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1 Introduction 

In an expanding universe distances grow following the expression 

r(t)=r a{t), (1) 

where a (t) is the scale factor, which indicates how a distance r varies with time, with ro being the present 
value of r. A Taylor's expansion of such expression around the present time to produces the approximate 
relation 

a (t) « 1 + H (t - t ) - -q Hl (t - t f , (2) 

where a (to) = 1, Ho = (a/a) \t=t and qo = — (da/a 2 ) \t=t > with a meaning the time derivative of 
the function a. There are several ways of measuring the quantities Hq (called Hubble constant) and qo (the 
decceleration parameter). Presently, the most cited method consists of the observation of supernovae in distant 
galaxies [1,2], and its use brought to light the idea that the universe is passing through a phase of accelerated 
expansion. So, a major problem confronting cosmologists today may be resumed in the question "what causes 
the acceleration of the universe's expansion?" Since there are yet no definitive and compelling answers, such 
issue, named as the problem of dark energy, remains open to debate. 

Another parallel puzzle of modern cosmology involves the contradiction between the observed movements 
of stars in the periphery of galaxies (and of galaxies in clusters of galaxies) and the movements expected from 
the amount of matter observed directly in such systems. This problem, also noticed when one analyses the 
data obtained from the cosmic microwave background 0> was christened as the problem of dark matter. 

Therefore, today's standard models of cosmology must face the above two problems and, therefore, the 
models used possess in general five basic components: photons, baryons, neutrinos, dark matter and dark 
energy. The question "what are dark matter and dark energy?" has, in such models, a standard answer J4): 
weakly interacting massive particles (WIMPs) combined either to a cosmological constant (yl-CDM models) 
or to some kind of scalar field (quintessence models). Both cosmological constant and a quintessential field 
could be the source of energy for the acceleration of the universe, while the weakly interacting massive 
particles, non-luminous, would be responsible for the effect of dark matter. 

Another kind of proposal consists in supposing that maybe dark matter and dark energy are in fact different 
manifestations of a single exotic entity, a component named quartessence (or unified dark matter, UDM) 0- 
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[32 1 . Quartessential models possess a single fluid whose main characteristic is to present an exotic equation 
of state p = p(p), which leads to a positive pressure in early phases of the universe and to negative pressures 
in the late phases. In such context, the main goal of this article is to present a study indicating how density 
perturbations can evolve in a universe composed solely of a single quartessential fluid. Here, the specific 
model used is known as the "modified Chaplygin gas" H331 - H361 . It is very important to notice that solutions 
for the scale factor obtained in this model may be obtained also in another kinds of cosmological models 1371 . 
Therefore, the study of the evolution of perturbations in each one of these models is an important step in the 
direction of 'breaking the degeneracy' between them. 

The structure of this work is the following: first, the modified Chapligyn gas model is presented, followed 
by another section with a brief - and somewhat pedagogical - review of the theory concerning the evolution 
of perturbations, leading to a differential equation which is used, in sequence, to show how the perturba- 
tions evolve in classical single fluid models. In the main section of the article the evolution of perturbations 
for a modified Chaplygin gas model [38] is obtained, with conclusions being presented in the last section. 
Throughout the text natural units are used, where c = G = % = ks = 1, unless stated otherwise. Also, 
arbitrary constants present in mathematical solutions are always labeled as c\ and c^. 

2 The modified Chaplygin gas model 

The modified Chaplygin gas is characterized by the equation of state 



where 7, M and p are free parameters. For 7 = // = 1 one has the original Chaplygin gas, while for M = 
one obtains the usual linear equation of state. One must notice that such gas is an ideal fluid, in the sense that 
it does not present variations in space, is, p — p (t). Therefore, the modified Chaplygin gas model does not 
incorporate any anisotropic stress, what distinguishes it from some more elaborate models [39|. 

Using the above equation, together with the condition for conservation of energy in an expanding universe, 



p= (7 -l)p- Mp~^ 



(3) 



p + 3- (p + p) = , 



(4) 



one obtains the relation 




(5) 
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which, when applied in the Friedmann equation, 



k 8tt 



-p, 



(6) 



yields general analytical solutions, valid for any value of the curvature parameter k, only for some values of 
7 and p. For example, for a flat space, Debnath, Banerjee and Chakraborty [ 34] give the general formula 

M 



Po^j-M \a 



\ 3 7 (1+m)' 



where a; = 1/ [2 (1 + p)], which can be inverted easily to yield a = a (t) only for /i = and p = —1/2, while 
for spaces with any curvature, explicit solutions for a (t) for the particular choices p = —1/2, 7 = 2/3 and 
7 = 4/3 are given by Costa J38). 

It is usual to choose p > 0, but here we will restrict the analysis to the particular choice p = —1/2, when 
one has explicitly, 

p=( 1 -\)p-Mp 1 l 2 , (8) 



what gives, for a flat space, where the curvature is null, 



+ [Po 



1/2 



M 



(ID 



3ff 2 



and, consequently, 



a(t) =X(e at -l) 



(9) 



(10) 



where A is a constant linked to 7, M, po and ao, and a = V 6irM 2 . Therefore, in such model 

2A 



_ a 2 e at f„at 



( 6 --l)*- 2 (^-l) 



(11) 



37 \37 
implying that if e Qt > 37/2 one will have a > 0. The conclusion is that values of 7 greater than 2/3 are more 
realistic in this model, since then the universe will experience an acceleration on its expansion only in later 
moments. 

Observing equations (HJl and (0! and considering the final argument of the previous paragraph, a 'good' 
choice could be 7 = 4/3, since in such case the equation of state would interpolate from a phase similar 
to a radiation dominated universe (p oc a~ 4 , for a << 1) to a phase with negative pressure, causing the 
acceleration of the expansion (for p < 9M 2 ), presenting also a phase of pressureless matter (since for values 
of p sa 9M 2 , one has p w 0). In resume, then, this particular modified Chaplygin gas model, with p = —1/2 
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and 7 = 4/3, could, at least in principle, play the role of an effective equation of state for a mixture of 
radiation, matter and a dynamical cosmological constant. 

It is very important to remark that the solution for the scale factor given by eq. ( TTOl) (the background 
solution) can be obtained also in another kinds of flat cosmological models, such as in a decaying vaccum 
cosmology II40II . Since an important difference between these two models is on the equation of state used 
in each one of them, a cosmological test based on properties affected directly by the equation of state could 
be a potential way to distinguish the models. The analysis of the evolution of perturbations, which has a 
dependence on the equation of state, is, therefore, a possible way to characterize a model in a way enough 
to separate it from its concurrents. In the next section, then, we present an overview of the general theory 
used for describing the evolution of perturbations of density, in order to show the particular relevance of the 
equation of state in such approach. 

3 The theory for the evolution of perturbations 

To make a comprehensive overview of the theory concerning the evolution of perturbations in an expanding 
universe, it is safe to begin with a Newtonian approach, where the radial acceleration R suffered by a test 
particle previously in equilibrium on the surface of a homogeneous sphere of radius R, due to the gravity 
caused by the increase in the sphere of an excess of mass Am, is given by the expression II41II 



where p is the mean density of the sphere and 5 is a density fluctuation (known also as density contrast), such 
that the density inside the sphere is given as p = p (1 + 5), with 5 « 1. The total mass of the sphere is 
m = Anp (1 + 8) /3, from where one can write, approximately, 



Using the relations given by equations ( 1121) and ( 1131 ), and noticing that R w R , one can write, then, 




(12) 




(13) 



R 
R 



6 



(14) 
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and, therefore, 



6 = AnGpS . 



(15) 



6 



This is the equation governing the evolution of density perturbations in a Newtonian static universe. However, 
if the universe is expanding, one can generalize the above arguments to obtain an equation where the scale 
factor appears explicitly 1411 . 

5 + 2-5 = AnGpS. (16) 

a 

The procedure leading to a fully general relativistic equation having the density contrast 5 as variable is 
a more elaborate one. Its development can be traced back to Hawking [42], with further advances by Olson 
[43 1 and, more recently, by Lyth and Mukherjee [44 1, and Lyth and Stewart [45 1. A resume of it can be found 
in a textbook by Padmanabhan H461 . Here we just cite the final result, which is, explicitly, for a Friedamnn- 
Lemaitre-Robertson-Walker (FLRW) fiat universqj (where the curvature parameter is zero), 



5 + H [2 - 3 (2w - v 2 )} 5 - ^H 2 (l - 6v 2 - 3w 2 + 8w) 5 = ~k 2 ^5 , (17) 



where H = a/ a is the Hubble parameter, w = p/p is the ratio between pressure and energy density of the 
cosmological fluid, v 2 — J 2 is the square of the sound velocity in the fluid, and k is the wavenumber of the 
Fourier mode of the perturbation. Notice that now one has a dependence on the equation of state p — p(p) of 
the cosmological fluid in study, included implicitly in the definitions of oj and v 2 . 

To better understand the processes described by these differential equations one can use a purely pragmatic 
point of view to see whether they can be obtained from effective "Lagrangian densities" where the degrees of 
freedom are given by 5 and 5. From the Lagrangian density 

C = g + 2,Gp5 2 ^ = \ + 2npc 2 5 2 (18) 

one produces the Newtonian equation (115b , with the constants c and G being used to give the correct di- 
mensionality. In this Lagrangian one has a usual kinetic term, proportional to 5 2 , and a negative potential 
proportional to S 2 , which implies there are no stable points of equilibrium. The solutions, thus, represent only 
decaying or growing modes. 

The result for an expanding universe, eq. ( 116b , can also be obtained by a Lagrangian densitjc 



C = C a a 2 = -\~— + ~-^-y (19) 

1 See Appendix [A] for another forms of this result. 

2 One way to build this result is to assume that the Lagrangian density given in equation dTsT > is an invariant quantity. Since the 
fraction c 2 /G is invariant, and distances grow in an expanding universe, while frequencies are shortened, the remaining quantity, 
which has dimensions of a squared frequency, must be multiplied by the square of the scale factor. 



c 



2 a 2 5 2 3a 2 H 2 5 2 
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with H 2 standing in the place of 8irGp/3. Again the potential is negative, but since one has now the kinetic 
term proportional to a 2 S 2 , there is a 'friction' in the system. 

It is not hard to verify that the general relativistic equation (Y7\ . valid for a flat space, can be obtained 
from a generalization of the above Lagrangians, 



where g is the determinant of the metric, and 1Z = K 2 — KijK' 1 ^ — 3 1Z is the Ricci scalar built from the extrinsic 
curvature tensor Kij and the three-curvature 3 1Z 1471 . For a flat FLRW model g = —a 6 and 1Z = 6H 2 . Again 
the potential is proportional to 5 2 but it can be positive or negative, according to the value of a polynomial 
expression involving u> and v 2 . Therefore, in the general relativistic case one cannot say that all solutions have 
only growing and decaying modes, since one can have oscillatory behavior. In order to make more clear this 
point, we present some examples in the next subsection. 

3.1 Perturbations in classical models of a single fluid 

As commented before, the evolution of perturbations in the general relativistic case can present three distinct 
possible behaviours: growing, decaying and oscillatory. This can be seen in the simplest models of cosmology, 
which use a linear equation of state, 



where 7 7^ is a numerical factor. Such relation, when applied into the energy conservation equation, 




(20) 



P= (7- !)/>, 



(21) 



p + 3H(p + p) = 0, 



(22) 



and in conjunction with the Friedmann equation for the flat space, 




(23) 



yields, for the scale factor, 




(24) 



and, consequently, 



2 1 

37 i ' 



(25) 
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Using these results into equation dl7l) . one then ends with the equation 

2 , v X 2 / n \ _ 



3- 



37t 3~fH 2 a 2 t - 

For any value of 7 7^ this equation has analytical solutions, since it is a particular case of equation 8.491.12 
from Gradshteyn and Ryzhik l48ll . 

z 2 u" + (2a - 2/3v + 1) zv! + [fi^z 20 + a{a- 2(3v)] u = , (27) 



with 



or 



2 2 fc 2 ^ . . / 2 x ~ 2 
3 7 a 2 V 3 7/ 

2 , /l 1 \ / 2 



4 , /l 1 \ / 2 x 1 



-2 H and v = — \ — I 1 

3 7 \2 3 7 y V 3 7 



3.1.1 Dust matter: 7=1 
In this case the equation to be solved is 

3t M 2 



$ + " = , (28) 



with the solution 

6 = ^r 1 + c 2 t 2 / 3 , (29) 

where the Cj are arbitrary constants. Here, therefore, the solution obtained contains only a growing and a 
decreasing mode. 

3.1.2 Radiation: 7 = 4/3 
In this case, one has 

with the oscillatory solution 



If 1 e k t0 

4t 



5 + —5 ^6 = %-5, (30) 

2t 2t 2 3 ah 



/Sin font 1 / 2 , 1/9 \ / COsbnt 1 / 2 .,1/91 

S = c 1 [ ° /2 - cos M 1/2 j + c 2 f ° /2 + sin M 1/2 ) , (31) 



where b\ = 4k 2 to (3a 2 ,) 1 and, again, the c, are arbitrary constants. 
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4 Evolution of perturbations in the modified Chaplygin gas model 

As seen in Section[2] the model we will study is defined by the equation of state ©, 

p={>y-l)p-Mp 1 '*, 

which has two free parameters, 7 and M. In order to verify how the perturbations in the matter density evolve 
in such model we need to calculate the quantities lu — p/p and v 2 = |^ which appear in equation ( 117b . 
Therefore, in the model analysed here one has explicitly 

2a 



w = V - = 7 - 1 - Mp- 1 ' 2 = 7 - 1 - 

p 3H 

2 d P 1 M -1/2 -, a 

V =dp- =1 - 1 -Y P ^^ , - 1 -3H 



(32) 
(33) 



Substituting the results above in the equation for perturbations, one has 

5 + [(5 - 3 7 ) H + 3a] 6 



1 22 \ 4 
(3 7 2 - 87 + 4) H 2 + ( — - 4 7 J aH + -a 2 

k2 (^n- 1^(7-1 



35 _ 
2 

Vis- 



(34) 



A 2 \2a 

It is easy to notice that this equation simplifies for 7 = 4/3. Then, this one will be the choice used in the next 
subsection. Anyway, one can rewrite the above equation using a new variable, y = 1 — = (1 — e at ) 1 , 



to produce 



. d 2 5 



10 \ 10 

4 } y 

37 J 37 



dS 



dy \l-y 



X 2 a 2 



i-vY 



1 



7 



67 2 - I67 + ; 
37 1 

1 1 



67 



3i 2 y 
7 1 



(5 = 



i-y y/ 2 (i-y) 



(35) 



For the particular case of the mode with k = 0, this last equation has a general solution for any value of 7 7^ 0, 

(36) 



1+— / \ — 

S( y ) = Cl V—^ 2 F 1 ( 2,1- — ;2- — ;« ) +c 2 ^- 

i-y V 3 7 ' 37' V 1-2/ 



where the c^s are arbitrary constants. 



4.1 Modified Chaplygin gas model with 7 = 4/3 



For the particular choice 7 = 4/3, equation (134b becomes 

<5 + (H + 3a) 5 - (H - 2a) (2H + a) 8 = - (l - J) (¥L - 1 ') ,) , 



(37) 
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In the particular case of the mode k = one has a solution given by eq. 

1 — p~ at f 1 "\ 1 \ 1 — e~ at 

i( " = c 'ITT^( 2 '2 : 2 ; — ) + «irr^f- (38) 

This solution, which is not oscillatory, can also be written in terms of the Hubble parameter H, as 

5 (if) = ckJi + c 2 5 2 , (39) 



where 



and 



«, - (40) 



(2H~af 1 /2# V /2iJ V 

^^S^ + hI-- 1 ) -tan(--ij , (41) 

and again the Cj are arbitrary constants. 
Approximate solutions 

An approximate solution valid for small values of t can be obtained directly from eq. d37b , using that a (i) « 

Aa 1 ^ 2 t 1 / 2 and H w t~ x ' 2 , so that the equation to be solved becomes 

<5 (5 fc 2 

5 + (1 + 6at) — + (4a 2 i 2 + Sat - l) — * = (l - 3ai + 2a 2 i 2 ) <5 , (42) 

it OA tt£ 

or, approximately, 

' + 5"5i— 53=*' (43) 

which is the equation for a classical model with 7 = 4/3, eq. (f30T >. with A 2 a taking the place of a 2 <g 1 . 

In order to obtain an approximate solution valid for large values of t it is interesting to rewrite equation 
03 as 



6 



«* - 1) 

Now, using that 5 (t) = e~ at (e at — 1) I 1 (i), the remaining equation is 

a , „s fc 2 



4 = -sk" (5^1 - 1 1 ' 1441 



3 v e- - zj J = - j^e (2 - e" 1 ) r , (45) 



which, for large values of at, can be approximed to 

A- 2 
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with the exact solutions 



r (t) = ci exp 



2^3 k 
3 \a 



+ c 2 exp 



2V3 k 
3 Xa 



valid for k ^ 0, and 



r (t) = ci — c 2 — e 



valid for = 0, where the c, are constants. 

Then, summarizing the approximate results for the modes with k ^ 0, one has, for t « 1, 

/sin fa 1 / 2 , 1/9 \ /cos fa 1 / 2 . , 1/2 

where now b 2 = 4k 2 (3A 2 a) , while, in the other extreme, for e at >> 1, 

/ 2^3 k _a± 3aA / 2^/3 k _ M 3ai X 

5t>>1=Cl exp^— _e"> -_j +C2 exp^- — _e" - - — 

with the Ci being arbitrary constants. 

4.7.2 Solutions in terms of scale factor and conformal time 



(47) 



(48) 



(49) 



(50) 



Although it is more interesting to see the density contrast as a function of the cosmological time t, it may be 
also useful to study its evolution in terms of other variables, such as the conformal time 77 or the scale factor 
a. The Appendix |Aj at the end of this work, shows such analysis for the classical single fluid models. Since 
in those models there is a substantial simplification of the differential equations for the density contrast, one 
can see whether or not this is also the case in the model of a modified Chaplygin gas. 

We choose to begin this extra analysis using the scale factor as the independent variable, when the relevant 
variable transformation is 



dt = 



2a 



a (a 2 + A 2 ) 



da , 



(51) 



and, therefore, the equation to be solved becomes 



, 2 ,0, d 2 S dS 4fc 2 a 2 — A 2 

(a 2 + A 2 ) —z + 8a— - —= 58 = -2 1 - ., 

v > da 2 da 3a 2 ( a 2 + A 2 ) V a 2 



A 2 



6a 2 + A 2 - a* 
a 2 + A 2 



5. 



Using the hypothesis 



A 2 

a (a 2 + A 2 ) 



r 



(52) 



(53) 
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one gets 



which for a — > yields 



12 ^2\ d r 2,3 , 2 * gLT 4fc 2 a 2 — A 2 _ 

(a + A ) — -5- + - (a 2 - A 2 ) — = — r *r, 

v ; da 2 a v ; rfa 3a 2 ( a 2 + A 2 ) 2 



d 2 r 2 dr 4fc 2 



da 2 a da 3a 2 A 4 



r = o. 



. sinfeia , \ / cos6ia . , . 
1 a ^o = c i — ; a cos Oia + C2 — : h a sin Oia , 



with the solution 

where b\ = 4k 2 / (3a 2 A 4 ), while for a >> A one has 

d 2 T 2 dr 4k 2 



da 2 ^ a da 3a 2 a' 



-r = o, 



with the solution 



A»A = c ie fclA2 / Q + c 2 e- blA2 / a . 



Finally, for fc = there is an exact solution, 



A=o = Cl 



1 /a\ 
— arctan — 
2A \\J 



2 (a 2 + A 2 ) 



+ c 2 . 



Therefore, in resume, one has 

c\ f sinfria 

Oa^O 



(a 2 + A 2 ) V ha 



cos 6ia 



C2 ( cos 6ia 



(a 2 + A 2 ) V ha 



sinbia 



>>A 



and 



4=0 = ci 



a (a 2 + A 2 ) 
A arctan (y) 



A 2 



A 2 



c 2 



a (a 2 + A 2 ) 



2a(a 2 + A 2 ) 2(a 2 + A 2 ) 2 
We can also present equations for evolution of the perturbations in terms of the conformal time, 

2 



V 



I ( -ai\ 

- — arccos e 2 
Aa V / 



so that 



a = A tan 

2 



(54) 



(55) 



(56) 



(57) 



(58) 



(59) 



(60) 
(61) 

(62) 



(63) 



(64) 



with < rj < f . With such substitutions equation (154b . for example, becomes 
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with the solution 

r - Cl 
r k=0 - - 

for k = 0, which is a result completely equivalent to equation d59l ). 

Equation ( 165 b produces a final important result, which can be obtained changing the coordinate to z = 
e~^- = cos {^^j • Doing this, the equation becomes 

( 1 - Z V + ^ = ^( 1 - 2Z ) f ' (67) 
which is a spheroidal wave equation of the type 

(1 _ z 2) ^ _ 2 ( a + 1) + (b - c 2 z 2 ) w = 0, (68) 
v ' dz z dz v ' 

as defined by Stratton H49II501 . with a = — 3/2, b = — and c 2 = —2b. Solutions to this equation are 
usually given in terms of elaborate series of special functions or orthogonal polynomials, such as 115 111521 

w(z) =]Td„T,?(z) , (69) 

n' 

where the e?„'s are the cofficients of the expansion, and T™ (z) are the Gegenbauer polynomials related to the 
hypergeometric series 2 Fi by 

™, / s (n+ 2m)! / 1 - z\ 

T n (*) = 2 m n \J. 2Fl r + 2m + h m + 15 ■ (70) 

The series d69l l is convergent only for a certain discrete set of values of the quantity A = b + a + a 2 , values 
which can be obtained from a transcendental equation 1531 . and which can be ordered so that for a given a, 
the lowest value of A is labeled A a . a , the next is A a <a +i, generating an index £ associated with the solutions 
through the eigenvalues A a j < A a £ +1 , and so on. The prime on the summation indicates inclusion of only 
even n's when the quantity I — a is even and only odd n's when £ — a is odd. Various normalization schemes 
can be used for the c? n 's, which satisfy a three-order recurrence relation 1541 . 

The subject of spheroidal wave functions is a particularly complex one, and therefore we will, in sequence, 
resort to numerical methods in order to better visualize the general behaviour of the density perturbations. 



sin (Xaij) 
Xa 



+ c 2 , 



(66) 
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4.1.3 Numerical solutions 

Besides the analytical solutions presented in this article, we also performed a numerical integration of Eq. (!3Vb 
in order to analyze the general behavior of the perturbations. The method used for the numerical calculation 
of the second order differential equation was the finite difference method. The code used for the numerical 
calculations was tested with the analytical solution obtained with k = 0, say Eq. d3 8 b - 

We perform the integration in time of Eq. ( f37T) starting from t — 10 2 years to t = 1.2 x 10 11 years. The 
approximate solution (|49| l was used to apply the initial conditions of our problem. The numerical results show 
clearly the presence of a peak in the evolution of the perturbations for any value of k. Results for some values 
of the wavenumber k are depicted in Fig. Q] The integral path used in Fig. Q] was At = 1 year. Numerical 
simulations with less acuracy, At > 1 year, show that the time in which the peak is maximum is practically the 
same but some details on the perturbation profiles are lost. Note that the maximum value of the perturbation 
increases when we increase the wavenumber value. The time in which the perturbation value is maximum 
also depends on the value of k. This behaviour is depicted in Fig. [2] The numerical results of Fig. [2] were 
obtained with a integration path of At = 10 3 years to save integration time. All figures were built with the 
values a — 1.28 x 10~ 10 years -1 , t = 1.4 x 10 10 years and a — a (to) = 1. 



k = 0.021 hMpc 
k = 0.022 h Mpc 



to 




5.0xlo'° l.Oxlo" 
t [years] 

Fig. 1 Graph showing the general behaviour of perturbations in a modified Chapligyn gas with 7 = 4/3 for two values of the 
wavenumber k. Here, it is clear the presence of peaks. The scale in the vertical axis is arbitrary. Numerical results for other 
wavenumber values present similar behaviours. 
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Fig. 2 Graphs presenting the dependence of the time in which the perturbation is maximal, t pea k, with the wavenumber of the 
perturbation, fc, for 7 = 4/3 (left) and 7 = 2/3 (right). 



4.2 Modified Chaplygin gas with 7 ^ 4/3 

Another values of 7 different of 4/3 allow one to obtain some analytical and numerical results, using proce- 
dures analogous to the one used for 7 = 4/3. Starting from equation ( 1341 ) or d35l l, using as variable the scale 
factor a, and putting 



6 = 



(-1)* i-af 



r. 



one obtains 



o\^l d 2 r 3 7 



da 2 2a 



dr 9-y 2 k 



2 1.2 



da 4a 2 a 4 



!f) 37 



37 

:*) t 



-1 2 



As a first example, for 7 = 2/3 one has then 



(71) 



!f) 37 



i + (x) 2 



r. 



(72) 



d 2 r 

(a + X)-^ + (a-X) 



1 dr 

a da 



k 2 X 



(a + XY X(a + X) 



r. 



(73) 



which for a — > reduces to 



d 2 r idr k 2 r 

da 2 a da 3a 2 A 2 a 2 



(74) 
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with the solution 



r a ^a = a 



c 2 a 



while for a » A one has 



d 2 r idr 2k 2 r _ 

da 2 a da 3a 2 a 4 



(75) 



(76) 



with the solution 



r a>>x = cxKq ( J +c 2 h 



(77) 



As another example, one can use 7 = 2 to obtain the equation 



(« 3 + a3 ) + - ( fl3 - A ') T" 

da z a 'da 



9k 2 X 



2\3 



(a 3 + A 3 )' 



(78) 



When a — > such equation becomes 



d 2 r 3 dr 9k 2 



da 2 a da a 2 A 6 



a 2 T = . 



(79) 



which has the solution 



n 2 



Ci Ji 



3fc 
2aX 3 ' 



a 2 +c 2 J- 



3fc 
27a 3 ' 



(80) 



while for a >> A one has 



d 2 T + 3dT + 9fc 2 A 3 r 
da 2 a da a 2 a 7 



= 0, 



(81) 



which has the solution 



r a »x = a 



ClJg 



6k\i 
1 

5a 



c 2 J_ 



6k\z 
5a 



(82) 



As numerical example, we choose the value 7 = 5/3 and, following the procedures presented for 7 = 4/3, 
we obtained the graphs presented in Fig. [3] which represent well the kinds of behaviours we encountered in 
our computations for several values of 7. 

We also made, now for 7 = 2/3, a numerical calculation of the time t pea k in which the maximum value 
of the peak is attained for each mode with wavenumber k, and such calculations are presented in Fig. [2] 
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Fig. 3 Graphs presenting the behaviour of the modulus of the density contrast, \5k\, with the time t (in years), for some values 
of the wavenumber of the perturbation, k, for 7 = 5/3. Notice the presence of oscillations in the early times, along with the 
formation of a peak. The vertical scale is arbitrary in all graphs. Also, in all graphs the integral path used was At = 10 years. 

5 Final remarks 

In this work, we present the mathematical analysis of the evolution of density perturbations in a flat cosmo- 
logical model with a modified Chaplygin gas acting as a single component, in order to verify whether this 
kind of analysis can be used to distinguish such class of models from another cosmological models which 
have the same solutions for the general evolution of the scale factor of the universe. For example, the solution 

a(t) = \(e at - 1) 1/2 , (83) 

shown here, can also be obtained in models with a constant bulk viscosity or with a decaying /1-term of the 
type A = aH ll3~7l . Our goal was to verify to what degree such degeneracy between models could be broken 
from the theory alone, without resorting to observational data. And our work shows that the analysis of the 
evolution of perturbations can be sufficient, since it produces results which can characterize uniquely a model 
such as the modified Chaplygin gas. As an example of this, we can point that a preliminary analysis of the 
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evolution of perturbations in the 4-decaying model for the mode k — 11551 indicates a result completely 
diferent from the one given by equation ( 1361) . 

It is also important to notice that we have shown it is possible to find a complete analytical solution for the 
evolution of the density contrast, involving spheroidal functions and the conformal time rj, in the particular 
model studied by us, i.e., a modified Chaplygin gas model described by the equation of state, 

p=( 1 -l)p-Mp-», (84) 

with 7 = 4/3 and [i = —1/2. It is very unlikely that other classes of models present the same kind of 
solution and, as a matter of fact, it is very unlikely also that other Chaplygin gas models, with values of 7 
and /1 different from the ones chosen by us, will lead to analytical solutions. However, it is certain that the 
analytical solutions we found, although limited to a specific choice of parameters in a model, can be useful 
to test numerical codes, to provide some insight for more generic numerical solutions, and by themselves are 
interesting for pedagogical purposes. 

Another relevant point to be noticed in our work is that our numerical analysis, together with the analytical 
approximate solutions obtained, shows that the density perturbation has peaks in its evolution. Furthermore, 
the time in which the peak is maximum depends on the value of the wavenumber k. Also, the perturbation 
magnitude is very sensitive to the value of k. These properties may characterize the modified Chaplygin model 
and help to distinguish it from other models. 

To finish, we must emphasize that a thorough analysis of the evolution of perturbations in a particular 
model must include the power spectrum, in order to allow a fair evaluation of the possible validity of the 
model. Here, our intention was to show that the results obtained for a specific choice of the parameters of the 
modified Chaplygin gas are enough to allow the distinction from another models, from a purely mathematical 
point of view, i.e., even without the construction of the power spectrum. Therefore, a comparison of the power 
spectrum predicted by this model and the one observed remains as a point to be attacked in future works. 
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A Other forms of the equation for perturbations 



An interesting alternative form of writing the differential equation which describes the evolution of the perturbations in the 
general relativistic context, eq. { flTt , is obtained using the scale factor a as the variable, i.e. (cf. equation 4. 122 from Padmanabhan 
(46)), 



da 2 2a da y ' 2a 2 v ' 

Such result can also be written as 

+ W [l-3(2c-« 2 )]|-H 2 (l-6^-3 W 2 + 8 W )f 

where r\ is the conformal time defined by the relation dt = adrj, and 

1 da 



1 2 2r 

-k v o , 



H 



a drj 



For classical models of a single fluid in a flat space, with 7 7^ 2/3, one has 



where 7/0 is a constant, and 



a (77) = ao [ — 



2 1 



37 - 2 ry 

so that the equation obeyed by the perturbations is 

d 2 5 fj- %\ 2dS , f~/-2 \ 25 



_ = -*' (7-1)5, 



(85) 



(86) 



(87) 



(88) 



(89) 



(90) 



dr/ 2 y-y - I J r) dt] \ T — f / 
from where it is easy to see the simplicity of three particular cases (7 = 1,4/3 and 2) and which has a general solution in terms 
of Bessel functions, 

1 97-10 / , \ 

,5 (77) = 37-2 j l37+2 [rjkJ^ - 1 . (91) 

1 2 37-2 \ / 

On the other hand, using a as variable one has the equation 



da 2 2a \ 7 3J da U ; V 3/ 2a 2 a 4 " 3 ^ 



■3^1(7-1)*, 



which also can be solved with Bessel functions, 



97-10 

S (a) = a 4 J, 1 37+2 
* 2 37-2 



3fe / 7-1 V 37^2 



37 — 2 \^67rpoa 



37 



(92) 



(93) 



and from where it is easy also to see the simplification which can be obtained for some particular values of 7. 



20 



References 

1. Perlmutter, S. et al.: Astrophys. J. 517, 565 (1999) 

2. Riess, A.G. et aL: Astronomical J. 116, 1009 (1998) 

3. Bennett, C.L. et al: Astrophys. J. Suppl. 148, 1 (2003) 

4. Turner, M.S.: The Sciences 40, 32 (2000) 

5. Kamenshchik, A., Moschella, U., Pasquier, V.: Phys. Lett. B 511, 265 (2001) 

6. Bento, M.C., Bertolami, O., Sen, A.A.: Phys. Rev. D 66, 043507 (2002) 

7. Bento, M.C., Bertolami, O., Sen, A.A.: Phys. Rev. D 70, 083519 (2004) 

8. Makler, M., Oliveira, S.Q., Waga, I.: Phys. Lett. B 555, 1 (2003) 

9. Ananda, K.N., Bruni, M.: Phys. Rev. D 74, 023523 (2006) 

10. Ananda, K.N., Bruni, M.: Phys. Rev. D 74, 023524 (2006) 

11. Diez-Tejedor, A., Feinstein, A.: Phys. Rev. D 74, 023530 (2006) 

12. Giannakis, D., Hu, W.: Phys. Rev. D 72, 063502 (2005) 

13. Perrota, E, Matarrese, S., Torki, M.: Phys. Rev. D 70, 121304 (2004) 

14. Gorini, V. et al: Phys. Rev. D 72, 103518 (2005) 

15. Carturan, D., Finelli, F.: Phys. Rev. D 68, 103501 (2003) 

16. Zhu, Z.H.: Astron. Astrophys. 423, 421 (2004) 

17. Sandvik, H.B. et al: Phys. Rev. D 69, 123524 (2004) 

18. Scherrer, R.J.: Phys. Rev. Letters 93, 01 1301 (2004) 

19. Dev, A., Jain D., Alcaniz, J.S.: Astron. Astrophys. 417, 847 (2004) 

20. Avelino, P.P. et al: Phys. Rev. D 69, 041301 (2004) 

21. Amendola, L. et al: Journal of Cosmology and Astroparticle Physics 7, 005 (2003) 

22. Beca, L.M.G. et al: Phys. Rev. D 67, 101301 (2003) 

23. Velasquez-Toribio, A.M.: Braz. J. Phys. 36, 772 (2006) 

24. Amendola, L. et al: Phys. Rev. D 74, 063524 (2006) 

25. Amendola, L., Waga, I., Finelli, F: Journal of Cosmology and Astroparticle Physics 11, 009 (2005) 

26. Colistete, R., Fabris, J.C.: Class. Quantum Grav. 22, 2813 (2005) 

27. Reis, R.R.R., Makler, M., Waga, I.: Class. Quantum Grav. 22, 1 191 (2005) 

28. Reis, R.R.R., Makler, M., Waga, I.: Class. Quantum Grav. 22, 353 (2005) 

29. Reis, R.R.R., Makler, M., Waga, I.: Phys. Rev. D 69, 101301 (2004) 

30. Avelino, P.P et al: Phys. Rev. D 69, 041301 (2004) 

31. Makler, M., de Oliveira, S.Q., Waga, I.: Phys. Rev. D 68, 123521 (2003) 

32. Reis, R.R.R. et al: Phys. Rev. D 68, 061302 (2003) 

33. Benaoum, H.B.: |hep-th/0205140| (2005) 

34. Debnath, U„ Banerjee, A., Chakraborty, S.: Class. Quantum Grav. 21, 5609 (2004) 



21 



35. Chimento, L.P.: Phys. Rev. D 69, 123517 (2004) 

36. Barreiro, T., Sen, A.A.: Phys. Rev. D 70, 124013 (2004) 

37. Costa, S.S., Makler, M.: |astro-ph/0702418| (2007) 

38. Costa, S.S.: Braz. J. Phys. 35, 1046 (2005) 

39. Koivisto, T., Mota, D.F.: Phys. Rev. D 73, 083502 (2006) 

40. Carneiro, S., Pigozzo, C, Barges, H.A., Alcaniz, J.S.: Phys. Rev. D 74, 023532 (2006) 

41. Ryden, B.: Introduction to cosmology, Addison- Wesley (2003) 

42. Hawking, S.W.: Astrophys. J. 145, 544 (1966) 

43. Olson, D.W.: Phys. Rev. D 14, 327 (1976) 

44. Lyth, D.H., Mukherjee, M.: Phys. Rev. D 38, 485 (1988) 

45. Lyth, D.H., Stewart, E.D.: Astrophys. J. 361, 343 (1990) 

46. Padmanabhan, T.: Structure formation in the universe, Cambridge Universe Press (1993) 

47. Kolb, E.W., Turner, M.S.: The early universe, Addison- Wesley (1993) 

48. Gradshteyn, I.S., Ryzhik, I.M.: Table of Integrals, Series and Products, Academic Press (2000) 

49. Stratton, J.A.: Proc. Nat. Acad. Sci. 21, 51 (1935) 

50. Stratton, J.A.: Proc. Nat. Acad. Sci. 21, 316 (1935) 

51. Morse, P.M., Feshbach, H.: Methods of Theoretical Physics, McGraw-Hill (1953) 

52. Leaver, E.W.: J. Math. Phys. 27, 1238 (1986) 

53. Abramowitz, M., Stegun, LA.: Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, 
Dover (1972) 

54. Weisstein, E.W: Prolate Spheroidal Function, MathWorld - A Wolfram Resource, http:// math- 
world. wolfram.com/ProlateSpheroidalWaveFunction. html 

55. Carneiro, S.: Modelo cosmologico com decaimento do vacuo - novos resultados e limites observacionais, talk presented at 
the VI Workshop Nova Ffsica no Espaco (2007). 



